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ABSTRACT 

Extrasolar planets and belts of debris orbiting post-main-sequence single stars 
may become unbound as the evolving star loses mass. In multiple star systems, 
the presence or co-evolution of the additional stars can significantly complicate 
the prospects for orbital excitation and escape. Here, we investigate the dynam- 
ical consequences of multi-phasic, nonlinear mass loss and establish a criterion 
for a system of any stellar multiplicity to retain a planet whose orbit surrounds 
all of the parent stars. For single stars which become white dwarfs, this crite- 
rion can be combined with the Chandrasekhar Limit to establish the maximum 
allowable mass loss rate for planet retention. We then apply the criterion to cir- 
cumbinary planets in evolving binary systems over the entire stellar mass phase 
space. Through about 10^ stellar evolutionary track realizations, we characterize 
planetary ejection prospects as a function of binary separation, stellar mass and 
metallicity. This investigation reveals that planets residing at just a few tens 
of AU from a central concentration of stars are susceptible to escape in a wide 
variety of multiple systems. Further, planets are significantly more susceptible 
to ejection from multiple star systems than from single star systems for a given 
system mass. For system masses greater than about 2Mq, multiple star systems 
represent the greater source of free-floating planets. 

Key words: planet-star interactions, planets and satellites: dynamical evolu- 
tion and stability, stars: mass-loss, stars: evolution, stars: AGB and post-AGB, 
(stars:) white dwarfs 



1 INTRODUCTION 

Roughly one tliird of all stars in the Galactic dis c are 
components of multiple star systems ()Ladal |2006| ) and 
the majority of multiple star syste ms are thought to be 
binary systems (|Duauennov fc May or 1991). Further, as 
of February 2012, several tens of extrasolar planets have 
now been detected or are suspected of existing in binary 
systems. Some planets have been repor ted in systems 
with ev en higher stellar multipl ici ties (e.g. Cochran et al.l 
19971: iRaghavan et al.l l2006l: [Guenther et all l2009l: 



Beuermann et all [2OI0I. I2OIII : iKuzuhara et~all I2OIII : 
Potter et all l201ll : lQran et all [2OIII . |2012| '1. The latter 



Desidera et al.l bOllh . iMugrauer fc Neuhauseij (|2009l 'l 



provide a helpful list of known exoplanets in multiple 
star systems, as of 2009. In binary systems, a planet may 
orbit one of the star s in what is sometimes called an 
S-typ e orbit (see, e.g. Lowrance et al. 20021: Bakos et all 
120061 : lEggenberger et all I2OO6I : ICorreia et all 120081 ^ 
Alternatively th e planet may orbit both stars in a P-type 
orbit (see, e.g. ISigurdsson et all l2003l : iLee et all l2009l : 

* E-mail:veras@ast. cam. ac.uk 



case describes circumbinary planets, whose existence has 
been bolstered by the recent tr ansit-based discoveries 
of Kepler-16 b (iDovle et al.l l201ll '). and Kepler-34b and 
Kepler-35b l|Welsh et al.l l2012l l. Therefore, understand- 
ing the dynamics of planets in multiple star systems and, 
in particular, binary systems is becoming increasingly 
relevant. This understanding includes how such planets 
form, and how they die. 

One potential avenue for planetary death is dynam- 
ical ejection as the s tar evolves beyon d the main se- 
quence and loses mass (jVeras et al .I2OIII . hereafter. Paper 



I). E v idence for free-floating planet s (in iLucas fc Roche 
2OOOI : IZapatero Osorio et al.l l2000l . I2OO2I : iBihain et al 



2009') provides observational motivat ion to inve s tigate 
this physical mechanism. In particular. ISumi et al.l(|201ll ) 
recently discovered 10 wide-orbit or free-floating bod- 
ies and calculated that l.SlJ g free floating planets ex- 
ist per main-sequence star. Therefore, more planets may 
travel between stars than orbit them, and this vast pop- 
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ulation of free-floaters cannot be explained by planet- 
planet scatteri ng immediately subseque nt to system for- 
mation alone jVeras fc Ravmondll2012l '). In Paper I, the 
authors analytically described the conditions which can 
lead to planetary ejection from a single star. They as- 
sumed isotropic mass loss and demonstrated that three 
key factors must be taken into account. These are i) the 
mass loss timescale, ii) the planetary orbital timescale 
and iii) the total mass of the system. That study was 
limited to consideration of a single phase of stellar evo- 
lution for mass loss which occurred in a linear manner. 

Here, we perform a three-tiered extension to that 
work by considering i) multiple stars, ii) multiple evolu- 
tionary phases and iii) nonlinear mass loss. Our primary 
application is the determination of the prospects for cir- 
cumbinary (P-type) planets to escape as the parent bi- 
nary evolves. Doing so provides us with the foundation 
to understand systems of higher stellar multiplicities, sys- 
tems which often include tight binaries. 

The physical evolution of a single star can, in princi- 
ple, be described by just the star's initial mass and metal- 
licity. However, the phase space of binary star evolution 
is significantly broader. In addition to each star's mass 
and metallicity, the separation and eccentricity of their 
mutual orbit (possibly including orbital variations due 
to massive planetary companions) are additional param- 
eters which must be taken into account. The metallic- 
ities of both stars are often assumed to be equivalent 
because the stars are assumed to have formed from the 
same molecular cloud and not to have subsequently ac- 
creted enough planetary material to signiflcantly alter 
their original [Fe/H] values. Our understanding of the 
physics of stars is not yet good enough to definitively 
link the above parameters with the wind velocity, the ac- 
cretion rate on to the companion and the amount of spin 
angular momentum transferred between the stars. There- 
fore, physical properties such as these represent a further 
broadening of the potentially explorable phase space. 

We focus our study by considering a single planet 
which is far from the parent binary (which may be tight 
or wide) and by analyzing the amount of mass loss from 
the system as a function of time. In the absence of mass 
loss, the approximately elliptical orbit of the planet will 
be negligibly perturbed by the potentially complex in- 
teractions between both members of the binary. In ad- 
dition, such complex details of this star-star interaction 
need only be modeled to the extent that they correctly 
yield the mass- loss rate from the entire system. If the bi- 
nary separation is large enough so that both stellar com- 
ponents evolve independently, then the planet's dynam- 
ical evolution reduces to the two-body case which was 
presented in Paper I, with the difference that the total 
system mass here is larger because of the presence of the 
binary companion. 

In Section 2, we consider theoretically how planetary 
ejection is affected by multiple phases of nonlinear mass 
loss for systems of any multiplicity and establish a crite- 
rion for retention we use throughout the rest of the paper. 
In Section 3, we numerically simulate binary star evolu- 
tion to identify which types of systems are susceptible to 
planet ejection and quantify where planets must reside in 
these systems in order to remain bound. We compare the 



binary and single star cases, and discuss the implications 
for higher multiplicities, other types of orbits, and exten- 
sions to this work in Section 4, and conclude in Section 
5. 



2 CHARACTERIZING MULTI-PHASIC 
NONLINEAR MASS LOSS 

This section provides the analytic background which 
motivates the numerical simulations in Section 3 and 
presents formulas which may be applied to other investi- 
gations. 

2.1 Phases, Regimes and Stages 

Here we define our nomenclature. We adopt the same 
meaning of phases which is commonly used in stellar evo- 
lution studies. Stages are treated as subsets of phases, and 
regimes refer the evolutionary environment of a planet 
due to stellar mass loss. 

2.1.1 Phase Identification 

Realistic stellar evolution occurs across several phases. 
A sequence of phases qualitatively characterizes a star's 
history. We adopt the definitions o f phase in t he SSE 
l|Hurlev et al.ll2000l ') and BSE (Hurlev et al.ll2002l ') stellar 
evolutionary codes: 

— Low mass (M < 0.7Mq) main-sequence star 

1 — High mass (Af > 0.7Mq) main-sequence star 

2 — Hertzsprung gap 

3 = First giant branch 

4 = Core helium burning 

5 — Early asymptotic giant branch 

6 — Thermally pulsing asymptotic giant branch 

7 — Naked helium star main sequence 

8 — Naked helium star Hertzsprung gap 

9 — Naked helium star giant branch 

10 — Helium white dwarf 

11 — Carbon/oxygen white dwarf 

12 — Oxygen/neon white dwarf 

13 — Neutron star 

14 = Black hole 

15 — Massless remnant. 

We denote phase numbers by k. Within each phase, the 
mass-loss rate can vary. Therefore, detailed modeling of 
a particular system necessitates fitting mass-loss rates 
with a piecewise nonlinear model. However, doing so for 
a broad exploration of stellar evolution phase space is not 
feasible. 

2.1.2 Regime Identification 

Instead, we seek to construct nonlinear mass loss profiles 
by a sequence of linear approximations. A key benefit 
of this approach is that we can utilize analytic results 
from Paper I and hence gain a better understanding of 
what conditions must exist for a planet to be significantly 
perturbed or to escape the system. Assume the stellar 
mass-loss rate is constant and equal to —a, where q > 0. 
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In order to characterize the tendency of a system to eject 
a planet, we use the dimensionless mass loss index, ^I* 
such that 



3' 



mass loss timescale a 



orbital timescale 



2-K \lMQyT- 



lAU 



IMr, 



(1) 



where a and n represent the planetary semimajor axis 
and mean motion, and /i = + Mp, where 

represents the masses of all of the stars in the system and 
Mp represents the mass of the planet. We use the term 
planet to loosely describe a bound body of any mass that 
is not a star and does not perturb the orbits of the stars. 

A planet is said to be evolving in one of two regimes 
depending on the value of 3'. In the adiabatic regime, 
when *!/ <C 1, a increases but the planet's eccentricity, e, 
remains constant. In the runaway regime, when '5^1, 
the planet's semi-major axis continues to increase and e 
can now vary (increase or decrease) across all possible 
values. Therefore, a planet may escape only if Vl* > 1. 
This bifurcation point is not exact. It is a weak function 
of both e and the planet's true anomaly, / (see Paper 
I). In fact, may instead be defined as aT j \i, where T 
represents the planet's orbital period. In this case, 3' does 
not contain the factor of l/(27r) that is present in Eq. ([l]). 
We will henceforth refer to this factor as k. 

If the planet is evolving in the runaway regime, the 
eccentricity evolution is a function of the eccentricity at 
the start of that regime and the planet's location along 
its orbit. For example, if the planet is close to pericentre 
on an already highly eccentric orbit, then the planet is 
ejected immediately. Alternatively, if the planet is close 
to apocentre on the same highly eccentric orbit, then the 
planet is never ejected. If the planet is on a circular orbit, 
regardless of /, it is ejected when the star has lost exactly 
half of its mass from its value at the start of that regime. 



2.1.3 Stage Identification 

Now suppose that within a given phase, the system un- 
dergoes A'^ consecutive stages of constant mass loss evo- 
lution with mass- loss rates oi a\,a2, ...,Q.i, ...,Q.n . Then, 
based on the results from Paper I, a planet experiencing 
i consecutive stages of mass loss in the adiabatic regime 
will have a semimajor axis and eccentricity at the end of 
stage i of 



ai = ao 



and 



e; = eo. 



(2) 



(3) 



where /3i = fii/fio represents the percent of the original 
system mass that is retained by the end of stage i. Be- 
cause eccentricity remains constant, a planet cannot be 
ejected in this regime (unless the semimajor axis evolu- 
tion carries the planet out of the system), even if almost 
all of one or both stars' mass is lost. The mass loss index 
at the end of stage i is: 



(4) 



where k = 1 or 27r depending on if the orbital timescale is 
defined with respect to the orbital period or the mean mo- 
tion. Equation Q indicates that escapability of a planet 
is independent of the details of the intermediate stages 
of mass loss provided that both and a; are known. 



2.2 Stage-based Retention Criterion 

Here, we describe the theory behind the critical system 
mass fraction lost for a given semimajor axis (Sections 
2.2.1 - 2.2.4) and the planet's critical semimajor axis for a 
given system mass loss prescription (Section 2.2.5). Com- 
putation of the critical semimajor axis will be the focus 
of subsequent sections. 



2.2.1 The Critical Mass Fraction 

Equation (|4]) demonstrates that, at the end of the ith 
stage of evolution, the minimum fraction of the original 
system mass which must be retained to guarantee that a 
planet with a given ao remains bound is 



/3crit 



= K3 



IMeyr- 



( — )' 

VIAU/ 



MO 

IMf- 



(5) 



Note importantly how this criterion can be written in a 
form which is independent of fii and a^. If /3crit; > 1, then 
the planet is not guaranteed to remain bound, regardless 
of how much mass is lost from the system. Equation ^ 
demonstrates that a single short but powerful ejection of 
mass from a star at any time can unhinge a planet. If 
this event occurs late in the life of a star, when is low, 
the probability is increased. Equation ((5} provides a use- 
ful way to quickly characterize an ensemble of systems. 
The equation holds for any stellar multiplicity. For sin- 
gle stars which become white dwarfs, the Chandrasekhar 
Limit can be inserted into fii. We may then determine 
the maximum possible Oi that can protect a planet. 



2.2.2 VisuaUzing the Criterion 

We can visualize the limits on a planet's escapability 
through Fig. [TJ which can be applied to systems of any 
stellar multiplicity in which a planet's orbit surrounds 
all of the stars. The upper and lower panels respectively 
illustrate planet retention bounds when one character- 
izes the mass loss index with respect to the mean motion 
(corresponding to k = 2tt; Eq. [l| and the orbital period 
— 1). The contours are values of /Scrit^ , and are a func- 
tion of ai and fio for a planet at ao — 10'* AU. Colours 
denote regions separated by adjacent contour levels in 10 
per cent intervals in masses. The upper panel of the figure 
demonstrates, for example, that a single 2M0 star which 
loses 60 per cent of its mass by the end of the ith stage 
of evolution is guaranteed to retain an orbiting planet at 
ao = 10'' AU if m < IQ-'^-^Mq yr-\ 
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Contours of yScrit, for ao = 10* AU and K=27t 




Contours of ;Scrit, for ao = 10 AU and k=1 



S -8 




Figure 1. Criteria to ensure planet retention. Each contour 
represents /^crit; i the minimum percent of the original system 
mass (^o) which must be retained by the end of an arbitrary 
evolutionary stage i in order to guarantee that a planet (at 
ao = 10* AU) remains bound. Different colours distinguish 
regions between 10 per cent differences in contour lines. In 
the upper panel, the mass loss index is defined with respect 
to the mean motion (corresponding to k = 2it; Eq. [T}. In the 
lower panel, the mass loss index is defined with respect to 
the orbital period (k = 1). The y-axis, which has the same 
scale on both panels, represents the mass loss rate from the 
system during the ith stage of stellar evolution. The thick 
black curve applies only to single stars where the planetary 
mass is negligible compared to the stellar mass and represents 
the maximum fraction of mass retained by the star if it is 
to become a white dwarf. Decreasing or increasing ao by one 
order of magnitude adds or subtracts approximately 1.5 units 
to the values on the vertical axes but otherwise leave the plots 
unchanged. 



Stars with masses up to 8Mq may become white 
dwarfs, whose maximum mass cannot exceed IAMq. 
Therefore, in the single star limit, the thick black curve 
denotes the Chandrasekhar Limit. Thus, in the upper 
panel, a 4Mq star which becomes a white dwarf has re- 
tained at most about 35 per cent of that mass during the 
transition. Planetary retention in such a system is guar- 
anteed only if the mass loss rate during the transition is 
^ IO-'^-^Mq yr~\ 

For planetary semimajor axes other than 10* AU, 
the plots are shifted vertically but otherwise remain unal- 
tered. For every order of magnitude that ao is decreased, 
the vertical axis values are increased by 10^'''' ~ 1.5 units. 
Therefore, for ao — 10^ AU, which represents a typical 
Oort cloud distance, Ui is rarely low enough to ever guar- 
antee protection of Oort cloud comets. 



2.2.3 Stellar Multiplicity Comparison 

We can also extend the figure to selected multiple star 
situations. Suppose that the primary has just become a 
white dwarf but the secondary is still on the main se- 
quence and is losing a negligible amount of mass such 
that M2 = = 'jMig. In this case, the maximum 

value of Oi which can guarantee planetary protection be- 
comes 

1 3 

1.4M0 -t-7Mio 



VaoMio (1 + 7) 



(6) 



This equation describes a limiting curve which can be 
drawn on Fig. [T]for a given binary system. Here, ai takes 
on an absolute maximum when jMig = 1.4Mq. Hence, 
the maximum possible value of the ratio of the binary to 
single star values of ai for a given Mig and 7 is about 
107 for the extreme bounds of Mig = 8Mq and 7 = 1. 
For the vast majority of realistic systems, however, this 
ratio is about a few. This exercise, which could be gen- 
eralized to higher multiplicities, demonstrates that the 
mass-loss rate for a main-sequence-white dwarf binary 
system needs to be slightly higher than in the single star- 
white dwarf case in order to eject a wide-orbit planet. 



2.2.4 Stellar Stage Comparison 

In order to deal with generic systems at all evolutionary 
stages, reconsider Eq. ((5|. Because /3criti is not explicitly 
dependent on the mass loss rates during the earlier stages 
of evolution, we can establish a critical value for the entire 
system evolution of 



^Scrit = Max (^c 



(7) 



where i — 1...N represents the ith stage of evolution. 

Comparing Pi with /3crit determines the prospects for 
ejection. If Pi > Pait for all Pi, then ejection cannot oc- 
cur. If Pi < /3crit for at least one pi, then ejection may 
occur. This is a general condition which may be applied 
to systems with any stellar multiplicity. In this work, we 
focus on binary systems, as they often represent compo- 
nents of higher multiplicity stellar systems and feature 
a representative set of physical processes that may be 
found in those systems. 
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2.2.5 The Critical Semimajor Axis 

Now suppose that is not given. Then we can deter- 
mine the critical value of ao, flcrit, for which the planet is 
guaranteed to remain bound through Eq. We obtain 



Ocrit 

lAU 



MO 



Merit ,z 



IMe/yr 



(8) 



where the minimum is taken over all values of i. Although 
IVeras fc WvattI (|2012l ) apply this criterion for the specific 
case of the Sun, the criterion is generally applicable to 
multiple star systems. We use Eq. ((8| to determine Ocrit 
throughout the remainder of the paper. 



2.3 Treating Nonlinear Mass Loss 

Ideally, a nonlinear mass loss profile may be approxi- 
mated with the finest possib le partition into linear seg- 
ments. IVeras fc Wvat^ (|2012l ) achieved this accuracy by 
considering the duration between every SSE timestep as 
a segment. Although effective for a detailed study of one 
system, this approach might not provide a consistent 
measure or be computationally feasible for a population 
of systems. Hence, we now consider alternatives. 



2.3.1 The Theory 

We compute Qfe by dividing the change in mass with the 
change in time over the entire phase fc. However, in cases 
where the mass-loss rate changes nonlinearly within a 
phase, we can model the mass loss evolution within the 
phase by consecutive stages of approximately linear mass 
loss. We select the boundaries for these stages based on a 
coefficient of variation, C„, of the data. If we fit a linear 
model to a given time series of data, then Cv represents 
the residual sum of squares about that fit divided by the 
mean of the response variable (the mass). 

An exact fit to the linear model yields Cv ~ 0. If 
we insist that each evolutionary stage satisfies Cv ^ C^, 
where Cv is a user-provided parameter, then the quality 
of our analytical treatment from Section 2 relies on the 
value of Cy. If Cv is too low, and a phase is split into 
too many stages, the resulting analysis becomes compu- 
tationally expensive. If Cv is too high, and each phase 
is represented by a single stage, then the analysis may 
represent a poor physical approximation. 

Suppose a phase fc is split into I stages, where j — 
Then Max(Qj) ^ because mass loss from the 
system is monotonic and at least one of the segments of 
a piecewise linear curve between both phase endpoints 
is steeper than the original curve. Therefore, the limit- 
ing value — > oo does not define any additional stages 
and hence provides a conservative estimate for planetary 
ejection. 



5.2M( 



0: 



2.3.2 An Example 

Consider the binary system with Mig 
Aho = 2.2M@, [Fe/H] = [Fe/H]o = 0.02, aso = WORe, 
and esf, = 0.0, where the subscript B refers to 
the primary-secondary binary. Suppose both stars 



began life on the main sequence. This system under- 
goes 9 distinct phases of evolution such that fci,2 ~ 
{(1, 1), (2, 1), (3, 1), (7, 1), (8, 1), (11, 1), (6, 1), (6, 15),' (11, 15)}. 
This sequence illustrates that Mi evolves off the main 
sequence first, on to the Hertzsprung gap and then the 
red giant branch. The primary's envelope is then blown 
away, leaving behind a naked helium star, which again 
begins a main sequence then Hertzsprung gap phase. The 
primary eventually becomes a white dwarf. However, at 
this point, the secondary is roughly three times more 
massive and is siphoning matter on to the primary. The 
primary is revitalized as a thermally pulsing AGB star, 
which soon after merges with the secondary. The new 
star continues AGB evolution until it becomes a white 
dwarf and dies out as such. 

The mass loss for the fci,2 ~ (6, 15) phase is markedly 
nonlinear. Fig. [2] shows the mass loss evolution of this 
phase for four different values of Cv. The latter three 
split the phase into stages. Phase and stage boundaries 
are given by vertical red dashed lines and brown dot- 
dashed lines, respectively. The total number of stages 
throughout the system evolution that are generated for 
Cv {oo, 0.1, 0.01, 0.001} are {9, 10, 12, 22}. The values 
of /Scriti are obtained from Ui, which in turn are com- 
puted from the intersections of the mass loss curve with 
the stage boundaries. Displayed in the lower-right panel 
are the three highest /3criti values of the many from that 
stage. Note that, as C* decreases, /3crit = Max (/Jcrit; ) in- 
creases. This must be the case because at least one of the 
segments of a piecewise linear curve between both phase 
endpoints is steeper than the original curve. The most 
accurate /3crit value is a factor of 1.86 greater than the 
least accurate value and corresponds to a factor of ~ 3.5 
in flQ. 



3 APPLICATION TO REALISTIC SYSTEMS 

In this section, we attempt to sample the entire phase 
space of binary systems in order to determine prospects 
for planetary ejection. We first argue that large regions of 
the phase space may be treated in a similar manner. For 
the remaining regions, we utilize a proven binary stellar 
evolution code that can generate tracks quickly, without 
solving coupled differential equations. Finally, we present 
our results through a series of contour plots that provide 
information about features of the binary systems studied 
as well as useful parameters such as the critical semimajor 
axis at which a planet is guaranteed to remain bound. 



3.1 Stellar Evolution Code 

We use a s lightly modified for m of the BSE stellar evolu- 
tion code (jHurlev et aklbOO^ ). This code utilizes empiri- 
cal, algebraic formulae derived from observational sources 
and theoretical frameworks in order to model evolving 
binary stars. Although the code tracks a wide range of 
stellar parameters, here we are concerned solely with the 
mass-loss rate from the system and the physical reasons 
for this mass loss. We assume that the circumbinary or- 
bit of the planet is always wider than the apocentre and 
common envelope radius of the binary. As shown in Paper 
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Figure 2. The representation of nonlinear mass loss as a sequence of consecutive linear stages. The coefficient of variation, Cv, 
about a linear fit to the data determines how to partition nonlinear mass loss within a phase. Plotted are partitions for the 8th 
stellar evolutionary phase of a system with Afi^ = 5.2Mq, M2g = 2.2Mq, [Fe/H] = [Fe/H]© = 0.02, a^^ = 100J?q and eg^ = 0.0, 
when Ml and M2 have merged and are undergoing AGB evolution. The thick gray curve represents the actual mass loss profile. 
The phase boundaries are given by the red dashed lines, and the stage boundaries within this phase are given by the brown 
dot-dashed lines. Values of /3criti are computed using Eq. l(5]l. The maximum of these values (Eq. [T]! determines the potential for 
the system to retain a planet and must increase as the number of stages increases. This potential is diminished as the profile is 
modeled more accurately. 



I, for monotonic stellar mass loss, a planet's semimajor 
axis must always be increasing, regardless of the shape of 
the orbit or the planet's position along that orbit. 

We have modified BSE with an updated form 
of the mass loss for stellar wind from naked he- 
lium stars. We have included a dependence on 
metallicity of M oc \fZ according to Eq. (22) of 
iNugis fc Lamerd (|2000l ). The mass-loss prescription 
we used is also incorporated in the NB0DY6 N- 
body code (which may be downloaded from the website, 
http:/ /www. ast.cam.ac.uk/~sverre/web/pages/nbody.htm 
as of December 2011). We also modified BSE to withhold 
output data dumps within evolutionary phases until 
at least 0.01 per cent of the entire system mass has 
been lost. This allows us both to use a mass-dependent 
uniform standard for analysis of all our simulations 
and to enable us to cover a wide region of phase space 
by limiting the number of outputs (particularly for 



main-sequence phases) in each simulation. We always 
output the first instance of a new evolutionary phase. 



3.2 Phase Transitions 

The phase and stage formalism described so far does not 
treat mass lost from a system during a phase transition. 
The BSE code computes the timescales and mass loss for 
most phase transitions. In these cases, we consider the 
mass loss between consecutive timesteps that encompass 
the phase transition to be linear. However, BSE treats 
some transitions as instantaneous. Although this approx- 
imation may be adequate for numerous stellar applica- 
tions, such as population studies, here that approxima- 
tion would cause a — >■ 00 and hence cannot be used. 

In the single star case, the only example of such an 
instantaneous transition occurs during a supernova. How- 
ever, for binary stars, there are several more examples. In 
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a comprehensive sampling of phase space, each of these 
instances must be accounted for and a timescale estab- 
hshed for each. We found that all such instances can be 
partitioned into four types of violent phenomena, i) core- 
collapse supernovae, ii) black hole core collapse, iii) ther- 
monuclear supernovae and iv) common envelope evolu- 
tion. In the following arguments, we assume that mass 
loss is isotropic and do not distinguish amongst neutri- 
nos, ejecta and gravitational waves as separate sources 
of energy or mass liberation. Accurate modeling of these 
phenomena would likely require removing these simplifi- 
cations and is beyond the scope of this study. 

3.2.1 Core Collapse Supernovae 

Commonly referred to as a Type II, Type Ic or Type Id 
supernovae, a core-collapse supernovae is a stellar explo- 
sion of a massive star (of 8 — 20Mq ) , typically while in a 
giant phase (fc G {3, 5, 6, 9}). The explosion leaves behind 
a neutron star (fc — 13) or a black hole (k = 14). This 
event may occur in a system of any stellar multiplicity 
and the fraction of the exploding star's original mass lost 
is typically 50 — 95 per cent. 

If we assume the tightest-known exoplanet orbit 
around a main sequence staiQEI (a ~ 0.015 AU) and the 
largest known progenitor mass which results in a neutron 
star (M ~ 2OM0), then no planet at any semimajor axis 
is guaranteed protection if a > O.9M0 hr~^. This con- 
servative rate is comparable to the mass loss timescales 
that are pr edicted from s upernova ejectio n velocities (e.g. 
iHamuv fc Pinto 2002; F esen et al.ll2007l '). 

Such planets are likely destroyed by the progenitor 
during its giant phase. For planets further away, with 
semimajor axes of 1 AU, this critical mass loss rate de- 
creases by a factor of about 544, and would be enhanced 
by just a factor of a few in a binary system, as argued in 
Section [2. 2. 3 1 Hence, realistic mass-loss rates from super- 
nova would well exceed the critical mass-loss rate for any 
planet that could survive the progenitor's pre-supernova 
evolution. 

Therefore, no planets orbiting stars of progenitor 
masses of between about 8 and 20Mq can be guaran- 
teed protection under the criterion of Eq. ([7]) and we can 
neglect sampling this region of stellar mass phase space 
in our simulations. There are some low-metallicity cases 
with 6 < M/Mq < 8 where the star might undergo core 
collapse supernova. In these cases, which correspond to 

fc„ G {0, 1, 2, 3, 4, 5, 6, 7, 8, 9} ^ {13, 14} (9) 

where ui = 1 or ui = 2, we assume that no planet is 
guaranteed protection. 

3.2.2 Black Hole Core Collapse 

Considering black hole formation in this work is par- 
ticularly important because all known black holes are 

^ See the Extrasolar Planets Encyclopaedia at 
| http: / /exoplanet .eu/ 1 

See the Exoplanet Data Explorer at jhttp: //exoplanets.org/| 



found in binary systems (iBelczvnski et al" ] |201ll ). Progen- 
itor stars with masses greater than about 2OM0 typically 
have extended layers of high density which prevent mass 
from escaping during an explosion. Nevertheless, the pro- 
genitor still forms a black hole through core collapse. The 
amount of mass loss which may accompany this core col- 
lapse is unconstrained by observations and can take on 
any value. A theoretical upper-bound for the t imescale 
for this process is 10 sec (|0'Connor fc Ottll20lil ). 

This mass-loss timescale is orders of magnitude 
shorter than the mass-loss timescale for core collapse 
supernova. Further, no planet at any semimajor axis is 
guaranteed protection in a system with the largest theo- 
retically postulated progenitor mass which might undergo 
core collapse and pr oduce a black hole (/i « 3OOM0, 
ICrowther et al.ll20lrf ) if Q > O.OI5M0/S. Therefore, any 
massive progenitor undergoing supernova- less core col- 
lapse need only lose a few tenths of a solar mass to place 
any orbiting planet in danger of ejection. In a binary sys- 
tem, this amount may be enhanced by a factor of just 
a few. Hence, we can claim conservatively that only in 
exceptional circumstances is a planet in this class of sys- 
tems guaranteed to remain bound. 

This claim allows us to neglect progenitor masses of 
M > 2OM0 in our simulations. 



3.2.3 Thermonuclear Supernovae 

Commonly referred to as a Type la supernova, a ther- 
monuclear supernova is an explosion and probable disin- 
tegration of a white dwarf (fc G {10, 11, 12} — ^ 15) ow- 
ing to mass accretion from a companion. Although this 
companion is often a giant star, the standard model has 
subgiant donors. Alternatively, the companion may rep- 
resent a white dwarf ( Pakmor et al..i2010i ). The velocity 
of the ejecta is comparable to that from core collapse 
supernovae but the to tal amount of mass ejected is typ- 
ically less than IMq (|Mazzah et al.l 120071 ). Despite this 
relatively low velocity of ejecta mass, the mass-loss rate 
is comparable to that from core collapse supernovae to 
within an order of magnitude. Therefore, because ther- 
monuclear supernovae afford little more protection for 
orbiting planets, we treat planets subjected to this type 
of mass loss in the same way as with core-collapse super- 
novae and assume that the planets cannot be guaranteed 
protection. 



3.2.4 Common Envelope Evolution 

Interactions between close binary stars resulting from 
Roche Lobe overflow may precipitate the formation and 
subsequent ejection of a common envel ope. Details o f 
common envelope evolution are complex l|lvanovall201ll ) 
and progress for greater understanding is being made 
thro ugh three-dimensiona l hydrodynamical simulations 
(e.g. lTaam fc Rickerllioiol : iPassv et al.ll201ll ). Therefore, 
establishing a timescale for common envelope evolution, 
tco, is difficult. However, this timescale, t^c, is thought 
to be within a couple orders of magnitude of the orbital 



© XXXX RAS, MNRAS 000.[TlfT9l 



8 Veras & Tout 



period of the binar^jf]. Even the lowest estimates for com- 
mon envelope evolution timescales guarantee that some 
planets remain bound. Therefore, this transition cannot 
be treated in the same manner as supernovae or black 
hole core collapse. 

In cases where BSE treats common envelope evolu- 
tion instantaneously, we assign a value to t^c- A value 
of tec = 10"^ yr represents an upper estimate for this 
timescale and with this value we can establish a con- 
servative estimate for escape. We also sample values of 
^ce = 10, 100 and 10* yr for relevant regions of phase 
space. 

3.3 The Lower Mass Extreme 

We have argued that planetary retention cannot be guar- 
anteed if at least one progenitor star has M > 8Mq . Now, 
we consider the lower mass extreme. 

The lowest-mass stars (M < 0.8Mq) have a life- 
time which exceeds t he current age of the Universe 
l|Parravano et al] [loij) . Their mass-loss rates are con- 
stant and negligible. By themselves, their mass loss evo- 
lution does not perturb a planetary orbit. However, in 
a binary system, their interaction with a more massive 
stellar companion could become a crucial aspect of their 
mutual evolution. Therefore, we include masses down to 
(M = O.2M0) in our simulations below. For stars with 
masses M <C O.IM©, planetary masses may no longer 
be comparatively negligible to the stellar masses. In that 
case, a less approximate treatment of the full three-body 
problem with mass loss would be required. 

3.4 Simulation Parameter Choices 

Having restricted the mass phase space to 0.2Mq ^ 
Mi,M2 S.OMq, reduced the nonlinear fitting of mass 
loss to one parameter, C^, and established the vari- 
able tec, we can now consider the binary's initial or- 
bit in more detail. We treat both close and wide bina- 
ries in this study, with binary separations of ub/Rq = 
10, 50, 100, 500, 10^ 5 X 10^ lO"* and 5 x 10'', and binary 
eccentricities cb = 0.0, 0.5 and 0.9. Wide binaries which 
do not interact with each other typically have separa- 
tions greater than 1O*-R0 « 47 AU). Therefore, in the 
extreme case of a_B = 5 x W^Rq, each component of the 
binary evolves off of the main sequence independently. In 
all cases, the planet is considered to be far enough away 
from the binary both to not affect the orbit of the binary 
and to have an approximately elliptical orbit around the 
binary. Each stellar evolution track is simulated for 15 
Gyr, and we adopt k = 2n in all our simulations. 

BSE contains several physical stellar parameters 
which may be varied. These include the accretion rate 
on to the secondary and mass-loss prescriptions. In order 
to focus our study, we use the default values for all phys- 
ical stellar parameters except for metallicity. For most of 
our simulations, we adopt either [Fe/H] = [Fe/H]© = 0.02 
or [Fe/H] — 0.0001, although we also sample metallicity 

^ Tills orbital period is commonly referred to as a dynamical 
timescale. 



values corresponding to [Fe/H] = 0.0005, 0.001, 0.005 and 
0.01. 

Throughout each stellar evolution realization, the bi- 
nary's orbit evolves. The components may interact with 
each other in complex ways, including merging through 
coalescence or collision, transferring mass in dynamical, 
nuclear and thermal regimes, rejuvenating a compan- 
ion through this mass transfer, and inciting accretion- 
induced collapse to achieve, for example, a supernova- 
less fc = 12 — > 13 transition. In all cases, the companions 
remain bound or merge. Only a supernova can unhinge 
the binary orbits. Common envelope ejection occurs be- 
yond the binary's orbit, and therefore should not affect 
the binary orbit. 

3.5 Simulation Results 

We present our results through a series of descriptive con- 
tour plots. On each, the x- and y-axes represent the initial 
mass of the primary and secondary star, respectively. We 
sampled each initial stellar mass at intervals of 0.2Mq, 
such that Ml ^ M2. Therefore, each plot represents a 
grid of 800 points. We can characterize many aspects of 
the system evolution with these plots, such as evolution- 
ary complexity, timescales and mass loss factors. Our pri- 
mary interest, however, is to determine the initial critical 
semimajor axis, Scrit, at which a planet is guaranteed to 
remain bound to the system. 

3. 5. 1 Dependence on Binary Separation 

The initial binary separation, aBo largely determines the 
sometimes complex manner in which binary components 
interact with each other when evolving off of the main 
sequence. Because the planet is considered to be far from 
the binary, if flcrit ^ as at any point, then the planet is 
not guaranteed to remain bound. The actual evolution of 
the planet in this case is complex and would necessitate 
detailed modeling. 

Figures[3]and[4]demonstrate how flcrit changes as asp 
increases from IQRq to 5 x lO^i?0 in both the [Fe/H] = 
[Fe/H]0 = 0.02 and [Fe/H]= 0.0001 cases. The dearth 
of black and blue contours on the upper panels of these 
plots immediately demonstrates how prone planets are to 
escaping these systems. 

First consider the black contours: all orbiting mate- 
rial is protected from ejection for M\,M2 ^ O.SMq in 
all cases, and that this protected region is extended to 
I.OM0 for the Solar metallicity case. For these lowest- 
mass stars, their proximity to each other is unimportant. 
The lower panels of the figures demonstrate that in this 
regime, stellar evolution undergoes only a few phases; in- 
deed, the stars never leave the main sequence. Even after 
15 Gyr, the stars do not lose enough mass on the main 
sequence to have a noticeable effect on orbiting mate- 
rial. This is the only evolutionary pathway that provides 
safety for Oort cloud comets with a > 10^ AU, which are 
predic ted to exist at galac tocentric distances beyond the 
Sun's (|Brasser et al.ll2010l '). 

Next consider the blue contours, found primarily 
in the widest initial binary separations, with ubq = 
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Figure 3. Fiducial prospects for planetary retention. The contours on the upper two panels represent the critical semimajor 
axis at which a planet is guaranteed to remain bound amidst binary stellar evolution as the initial binary separation is increased 
from left to right. Initially the binary has a circular orbit and solar metallicity, and we adopt tec = 10^ yr and Cv — > oo. The 
lower panel provides qualitative detail about the evolution of the systems with as = 50 Rq. The contour labels in the lower- 
loft plot are CE Only- at least one common envelope phase but no supernovae, SN la Only- thermonuclear supernova only, SN 
la+CE- thermonuclear supernova plus at least one common envelope phase, SN CoreCo Only- core collapse supernovae only, SN 
CoreCo-|-CE- core collapse supernovae plus at least one common envelope phase. Both SN- both core collapse and thermonuclear 
supernovae. The figure demonstrates that circumbinary planets are highly susceptible to ejection. 
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Figure 4. Low metallicity prospects for planetary retention. The upper two panels are the same as Fig. |3]but for [Fe/H] = 0.0001. 
The lower panel provides more detail on the the systems with aB = 500-Rq . The total fraction of system mass loss over its lifetime 
is a poor indicator of a^rit , whereas the total number of phases experienced by the binary hints at the reason for the vertical green 
strip at 3Mq ^ Mig ^ AMq. Although low- metallicity progenitor masses are prone to core collapse supernova for M > 7Mq, 
they afford circumbinary planets slightly more protection than higher-metallicity systems for wide binary separations. 
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5000i?Q « 23 AU. Comparing these two plots in Figs 
13] and 3] demonstrates tiiat the higher metaUicity case 
generaUy aUows for greater planetary protection unless 
■^lo ~ 6.6M0, when a core-collapse supernova occurs. 
A circumbinary planet with a < 1000 AU is guaranteed 
to remain bound if Mig, h'ha ^ 2Mq, although in many 
cases the planet may still be safe further away. This is 
the expected result from the single star case where mass 
loss is treated linearly in each phase (see Paper I), be- 
cause at this separation the stars effectively evolve in- 
dependently from each other. However, this treatment 
is conservative; later we demonstrate how one of these 
plots changes when CI is decreased and higher accuracy 
is achieved. For = IQQQ-Rq ~ 4.7 AU, planets with 
a > 1000 AU may still survive if Mi^ , ^ 2Mq , as 
long as Cu — >■ 00. The contour plots for hbq = 1 x lO^i?0 
and aBg = 5 x IQ'^Rq axe nearly indistinguishable from 
the aBo = 5OOO-R0 case and are not shown. 

For closer separations, the stars are more prone to 
interact with each other violently. The lowest-right plots 
in both figures hint at the complexity of t he stellar evo- 
lution for asa = 50 Rq and = 500RrT, . fHurlev et al.l 
(|2002l ) describe the complete details of the physical mech- 
anisms behind the variations in these contour plots. Here, 
we just point out some of the most important features. 

When ubq ^ IOOO-Rq, prospects for planetary pro- 
tection are primarily determined by common envelope 
and supernovae events. The lowest-left plot of Fig.[3]char- 
acterizes where in stellar mass phase space for ubq = 
5O-R0 these events occur. They fill almost the entire phase 
space. Any contours on this plot that indicate a super- 
nova has occurred automatically yield Ocrit < 5 AU (no 
protection) in the corresponding upper panel figure, as 
expected. Some contours which indicate that no super- 
nova occurred also yield Ocrit < 5 AU. In these cases, 
BSE computes a timescale for all phase transitions, and 
the mass-loss rate at some point during the evolution is 
great enough to prevent planetary protection. 

Now consider cases where 5 AU < Ocrit < 300 AU 
for ubq ~ 50-Rq. Because we adopted tec = 1000 yr for 
all these simulations, the variance in Ocrit is due to the 
amount of mass lost during the common envelope phase. 
Later, we will explore how the results change when tec is 
altered, especially to less conservative values. The green 
patch found at {4Mq < Mi„ < 8Mq, M20 < 2Mq} 
indicates that the fraction of system mass lost owing to 
common envelope evolution here is lowest. 

A conspicuous feature of the {aso ~ IOOO-Rq, [Fe/H] 
= 0.02} and {ubo = 500-Rq, [Fe/H] = 0.0001} contour 
plots is a vertical green patch surrounded by yellow and 
orange contours. This green patch arises because the com- 
mon envelope in these systems is formed at a different 
phase of stellar evolution than the surrounding systems, 
at a phase where less of the system mass is lost. 

Finally, the lower-left plot in Fig. |4]displays contours 
of the total fraction of system mass lost after 15 Gyr. 
The plot demonstrates that there is almost no correla- 
tion, except in the lowest progenitor mass cases, to the 
corresponding middle panel Ocrit plot. This comparison 
illustrates how the total system mass lost in a system 
is a poor proxy for Ocrit (also compare with Fig. 3 of 
IVeras fc Wvat?l2012l ). 



3.5.2 Dependence on Binary Eccentricity 

Related to initial binary separation is initial binary ec- 
centricity. The pericentre of the orbit helps determine 
how close to each other the stars orbit. Fig. [5] explores 
how the fflso = IOO-Rq and aso = IOOO-Rq Ocrit contour 
plots from Fig. [3] change when the binary's initial eccen- 
tricity is moderate [cbq = 0.5) or high [cbq = 0.9). The 
plots reveal that acrit is relatively insensitive to nonzero 
binary eccentricities unless the eccentricity is high. In the 
high-eccentricity case, the plots resemble those obtained 
from smaller initial binary separations an d circular or- 
bits. This is consistent with the finding of iHurlev et al.l 
(200 j) that it is the initial semi-latus rectum, or angu- 
lar orbital momentum, rather than semimajor axis and 
eccentricity which determine how a system evolves. 

However, there are some differences. For obq = 
100i?Q, with the initial pericentre at 10-Rq, a bifurca- 
tion in the plot occurs where Mip -I- M2q ~ 8Mq, re- 
sembling the circular IQRq case. However, the behaviour 
at {2Mq s; Mio ^ AMq , IMq sS Ma^ < ZMq), is 
markedly different. The red patch bounded by this re- 
gion in the circular case is caused by thermonuclear su- 
pernovae. However, in the highly eccentric case, the con- 
ditions favorable to producing this phenomena vanish. 
The resulting value of Ocrit is then dominated by com- 
mon envelope evolution. 

For osg = 1OOO7?0, with an initial pericenter at 
IOOJJq, the flcrit distribution is lower than in the circu- 
lar obq = IOO-Rq case. The explanation again relates to 
thermonuclear supernova. However, here the high eccen- 
tricity of the binary promotes the conditions necessary for 
this explosion, thereby suggesting that for high values of 
e_BQ, the dependence of Ocrit on bbq may be complex and 
best treated on an individual basis. 

3.5.3 Dependence on MetaUicity 

In Fig. |6l we explore how Ocrit changes during the tran- 
sition between [Fe/H] = 0.02 (Fig. [Sjl and [Fe/H] = 
0.0001 (Fig. mi through the ob,, = lOOOi?© case. As 
[Fe/H] decreases, primaries with Miq > 7Mq are more 
likely to experience core collapse supernovae. However, 
for Afio < 7M0 and [Fe/H] ^ 0.005, the plots are nearly 
identical. The primary differences occur for 4.0Mq < 
M\o ^ 7.OAf0 and h/h^ Q.GMq. In this regime, no 
supernovae occur, and Ocrit is dominated by common en- 
velope evolution. The amount of mass loss during this 
evolution is greatest for the higher metaUicity stars. This 
is the reason for the red strip on the bottom of the upper- 
right plot in the figure. 

3.5.4 Dependence on Common Envelope Evolution 

In order to demonstrate the dependence of common en- 
velope timescale on Ocrit, we select a set of systems for 
study whose stellar evolution is dominated by a common 
envelope phase that BSE treats as instantaneous. The 
upper left plot in Fig. [7] presents our selection. The other 
plots in the figure demonstrate how ttcrit changes when 
tec is varied from 10 yr to 10"* yr (the 10^ yr case is shown 
in Fig.[3ll. 
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Figure 5. Critical semimajor axis for planet retention as a function of initial binary eccentricity, sbq - The upper panel and lower 
panel can, respectively, be directly compared with the = lOOR© and as = lOOO-R© systems from Fig. |3] Moderate initial 
binary eccentricities (e = 0.5) have little effect on the systems. However, high eccentricities (e = 0.9) can cause close encounters 
which lead to phenomena that are associated with smaller values of as. 



Fig. [7] demonstrates how conservative our fiducial 
choice of tco = 10^ yr is. For most of the stellar mass 
phase space on these plots, the initial orbital period of 
the binary is less than 10 yr. Therefore, if common enve- 
lope formation and destruction proceeds on the order of 
a few dynamical timescales, few planets would be guar- 
anteed to survive. Those that do would have flcrit ^ isq' 
which destroys our assumptions and likely would destroy 



the planet as well. Additionally, note the sharp transi- 
tion between the blue and red contours in the tee = 10 yr 
case, indicating the importance of the onset of common 
envelope evolution. 
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Figure 6. Critical semimajor axis for planet retention as a function of stellar metallicity, [Fc/H] . These plots represent intermediate 
mctallicity values from the two extremes in Figs l3l and l4l for as = 1000-Rq. The critical semimajor axis is largely insensitive to 
[Fe/H] values below about 0.001. 



3.5.5 Dependence on Nonlinear Mass Loss 



Here we relax the conservative assumption C* — >■ oo for a 
subset of cases. Doing so for all stellar evolution realiza- 
tions would be ideal but is too computationally expen- 
sive. Further, because Ocrit is determined by phase tran- 
sitions for the vast majority of phase space, relaxing the 
assumption — >■ oo is largely unnecessary. However, for 



low progenitor masses and large separations, the choice 
of C* affects acrifl 

In Fig. [S] we choose a set of systems where flcrit 
is dominated by (upper-left plot). This set features 



* Unlike for multiple star systems, in single star systems with 
M < 8Mq , Ocrit is dominated by the choice of C* for the vast 
majority of the stellar mass phase space. 
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Figure 7. Critical semimajor axis for planet retention as a function of common envelope evolution timescale, tee- These plots 
encompass the representative range for the formation and ejection timescale of a common envelope (10^ yr in the upper right plot 
to 10^ yr in the lower right plot) and can be directly compared to the ag = IOOO-Rq systems in Fig. [S] for which the fiducial value 
of tee = 10^ yr is adopted. The upper loft plot demonstrates our motivation for choosing this set of systems for study: the common 
envelope phase dominates the stellar evolution mass phase space for the initial binary orbit. The contour labels in the upper-left 
plot are described in the Fig. |3] caption. 



widely separated stars which evolve physically nearly in- 
dependently of one another. We then computed the max- 
imum value of Cv achieved throughout each of the 800 
stellar evolution realizations (upper-right plot) and used 
those to sample values of Cy that would ensure splitting 
the phases with the most nonlinear stellar mass-loss rates 



into multiple linear stages for the vast majority of simu- 
lations. The — >■ oo case is plotted in Fig. O 

Setting Cy — 0.05 (lower-right plot) eliminates the 
dark blue contour from Fig. [3] for Mi^ > 2Mq and re- 
duces acrit by a factor of at least a few across the entire 
stellar mass phase space. This result reinforces the need 
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Figure 8. Critical semimajor axis for planet retention as a function of Cv The upper left plot, composed of 800 points, 
demonstrates our motivation for choosing a set of systems with ag = SOOOi?© : the vast majority of stellar evolution realizations 
featured in this set of systems feature well-defined and nonviolent phase transitions. Therefore, the critical semimajor axis is 
primarily determined by mass loss within phases rather than in between phases. The contour labels in the upper-left plot are 
described in the Fig. |3] caption. The upper-right plot demonstrates that our choices of C* = 0.10 and C* = 0.05 affect the vast 
majority of the systems sampled. The — > oo case is plotted in Fig. |3] Lower values of Cv provide slightly more accurate but 
significantly more computationally expensive results. 



for detailed modeling of individual systems that feature 
highly nonlinear mass loss during, for example, the last 
epochs of giant branch evolution. 



4 DISCUSSION 

Here we discuss four relevant related topics to this work. 
The first compares planetary ejection rates for single and 
binary stars, which helps inform their relative contribu- 
tions to the free-floating planet population. The second 
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importantly considers how the violent behavior which of- 
ten accompanies P-type planetary orbits extends to other 
configurations and higher multiplicities. The third dis- 
cusses the limitation of our knowledge of some aspects of 
stellar evolution. In the fourth, we consider partitions of 
the runaway regime into multiple stages, and show that 
in some cases there exist analogues to Eqs and @. 



4.1 Escape Rate Comparison With Single Stars 

Observations suggest that ne arly two free-fioat ing planets 
exist per main sequence star (|Sumi et aLlbOlll l . This vast 
population cannot be explain ed by instability-induce d 
planet-planet scattering alone (jVeras fc Ravmondll2012l ). 
Other potential sources of free-floating planets arise from 
escape in evolved single star and binary star systems. 
Although a detailed ejection rate computation based on 
initial mass functions, binary fractions, distributions of 
binary separations, and semimajor axis-based planetary 
distributions is beyond the scope of this study, we can 
provide some qualitative estimates here. 

The critical semimajor axis is given by Eq. which 
is applicable for any stellar multiplicity and a planetary 
orbit that surrounds and is far from the central star or 
stars. This equation contains the total initial system mass 
/io- On the critical semimajor axes contour plots in this 
paper (Figs. [3l8|l . lines of constant no would be approxi- 
mately diagonal from the lower right to upper left. There- 
fore, the tight binaries in Figs. [3][5] demonstrate that for 
Ho ^ 2Mq , a planet at several tens of AU away from the 
binary stars is typically prone to escape, and a planet at 
several hundred AU away is prone to escape in nearly all 
cases. In contrast, direct 2-body numerical simulations 
from Paper I illustrate that for single stars, for no value 
of ^0 does a planet with a semimajor axis under 100 AU 
escape. Further, at 500 AU, planets begin to escape at 
/io ^ 3Mq. Therefore, for a given system mass, planets 
are more prone to escape in a binary system than in a 
single star system. 

This conclusion is critically supported by the lack of 
thermonuclear supernovae (see Section 3.2.3) and com- 
mon envelopes (see Section 3.2.4) in single star systems. 
The bottom left contour plot of Fig. |3] illustrates that 
either or both of these violent phenomena occur in the 
vast majority of the stellar mass phase space for that 
set of tight {50Rq), Solar- metallicity binary systems. In 
particular, the common envelope which typically forms 
for 2 ^ fio/MQ ^ 6 values is not present in single star 
systems and promotes escape for all reasonable blow-off 
timescales (10 yr - 10'* yr). Even in moderately wider bi- 
naries, with binary separations of several AU, the stellar 
mass phase space is dominated by common envelope for- 
mation and ejection (see the upper left plot of Fig. [7| ■ If 
the common envelope blow-off timescale is as short as 10 
yr, then the bottom right plot of Fig.[7|demonstrates that 
nearly all planets, regardless of their semimajor axes, are 
prone to escape for /io ^ 4:Mq. In contrast, in the sin- 
gle star case, no planet under 100 AU will escape for 
Ho ^ 4:Mq (right panel of Fig. 14 in Paper I) . 

For the lowest mass systems, with total masses under 
one Solar mass, we can crudely estimate that planetary 



escape for single stars is negligible based on arguments 
from Paper I. Similarly, Figs. 1 3181 show that /io ^ ^Mq 
binary systems protect planets within 10^ AU. These 
low-mass systems represent the majority of stars in the 
Milky Way, based on a wide s electio n of initial mass func- 
tions from IParravaiw et all (120111 ). Therefore, we con- 
clude that these systems do not contribute to the free- 
floating planet population, regardless of m ultiplicity. For 
higher -mass systems, as summarized by iBastian et al.l 
l|2010l 'l. the binary frequency of Solar-type stars is ap- 
proximately 60%; this value increases to nearly 100% for 
0-type stars. Therefore, more stars appear in binaries 
than not for /io ^ 2Mq. Further, as argued above, plan- 
ets are more susceptible to escape from binaries than from 
single stars for /io ^ 2Mq. Therefore, assuming that the 
number of planets which survive main sequence evolution 
is approximately equal in systems of all stellar multiplic- 
ities, the contribution to the free-fioating planet popu- 
lation is greater in multiple star systems than in single 
star systems for /io ^ 2Mq. For the intermediate range 
1 Ho/Mq ^ 2, the comparison of ejection contribu- 
tions is less cl ear, because the bin ary fraction might be 
less than half l|Bastian et al.|[2010l ). 

4.2 Higher Multiplicities and S-Type Orbits 

A detailed extension of these results to planets on S-type 
orbits or to stellar systems with more than 2 stars is non- 
trivialfj. However, on qualitative grounds we claim that 
flcrit for both of those system types is generally lower than 
for the binary systems studied here. 

A planet orbiting one star in a close binary is likely to 
be destroyed before entering a regime where it is prone to 
escape. Although a planet may briefly survive engulfmen t 
by a stellar envelope of a single star (|Bear et al.ll201ll ). 
a planet residing in or near a binary common envelope 
would be subject to 3-body forces that would likely cause 
a collision. Similarly, a supernova by either star might 
incite a 3-body collision or escape if not direct destruction 
of the planet. A planet's orbit that is entangled in a mass 
transfer stream between the two stars would likely harbor 
a quite interesting but destructive orbit. 

If the planet orbits one star in a wide binary, then 
stellar evolution of the parent star would extend the 
planet's orbit, possibly beyon d a region of stability in 
the three-body problem (e.g. iHolman fc WiegertI 1 19991 : 
iDonnisonllioog ). If the parent star explodes, the planet 
is ejected, regardless of the mass of the binary compan- 
ion. If instead the binary companion explodes, then the 
planet and its parent star together will become unbound 
from the companion's remnant but could remain bound 
to one another. The planet's fate then depends on single 
star evolution. 

^ Kratter & Perets (2012) consider a planet wliicli can safely 
"hop" from one star to another during post-main-sequence 
evolution of an S-type binary and remain in a long-term sta- 
ble orbit after hopping. If, however, both of these stars are 
evolving on the main-sequence, then a planet may "bounce" 
between the st ars but fail to achieve a long-term stable orbit 
around either jMoeckel fc Verasll20l3) . 
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Now consider P-type orbits in systems with higher 
multiplicities. Examples include combinations of stars 
and possibly other planets which are all orbited by a 
distant planet in an approximately elliptical orbit. For- 
mation scenarios in which a planet can form from core 
accretion in a, for example, circumternary disc, have yet 
to be explored fully. However, second-gene ration planet s 
may form in a vari ety of exotic syst ems l|Peretsl [2OI0I ) 
and if the claim by ISumi et al.l (l201lh that free-floating 
bodies are more abundant than main-sequence stars is 
true, these planets may be captured by multiple star sys- 
tems. Although additional stars add to the total mass of 
a system, thereby inhibiting escape according to Eq. (O, 
this effect is overshadowed by the increased likelihood of 
repeated, violent mass loss events. Further, binary evolu- 
tion represents the simplest-possible outcomes for closely 
interacting stars, and violent phenomena that have not 
yet been characterized may exist for systems of higher 
multiplicity. Many multiple systems contain binary com- 
ponents and thus suggests that as a minimum an orbit- 
ing planet is subject to the restrictions suggested in this 
work. 



4.3 Improved Physics 

As researchers gain a better understanding of the under- 
lying physical processes featured by BSE, the quantita- 
tive results of this work are likely to be modified. Par- 
ticular events, such as common envelope evolution, are 
largely unexplored. Only recently has the dynamics of 
mass trans fer in binaries on non-circular orb its been in- 
vestigated. ISepinskv et~all (|2007l . |2009| . I2OIOI ) developed 
a formalism for mass transfer in eccentric binaries that 
represe nts a step forward in understanding this physical 
process. Ivan Rensbergen et al.l (|2008l . [2"oiOt ) find that the 
mass-transfer rate for low-mass binaries is never large 
enough to allow for mass loss from the system. They 
claim that matter can escape the binary if the kinetic 
energy from fast rotation plus the radiative energy of 
the hot spot exceeds the binding energy of the system. 
Mass transfer in cataclysmic variables might repeatedly 
vary by over an order of magnitude owing to nova out- 
bursts l|Kolb et al.ll200ll ). iMeng et all l|2008f ) considers 
the amount of mass lost for the highest metallicity stars 
([Fe/H] = 0.04 - 0.1), a regime not treated by BSE. Re- 
cent h ydrodynamic models of mass transfer between bi- 
naries (iLaioie fc Sillsll2011al lbh attempt to circumvent the 
limitations of the Roche lobe formalism, which was de- 
veloped for restricted cases. Although these simulations 
are just beginning to tackle the complexities of mass loss 
in eccentric binaries, they demonstrate that some trans- 
ferred material is ejected, while other becomes loosely 
bound. 



linear stages, we may derive equations analogous to Eqs. 
and (Ul. However, the evolution of a and e in the 
runaway regime is more complex (see Paper I) and we 
have derived analytical evolution equations only when 
the planet is assumed to begin and hence remain at peri- 
centre or apocentre. 

In the first case, for / — f{t) = 0° , despite the added 
complexity of the relations, they couple together to cancel 
out all intermediate stage terms so that ai and d may be 
expressed in terms of ao, eo and I3i only so that: 



gp (1 — ep) 
2 - (1 + ep) 



and 



=/?-^(l + eo)-l. 



(10) 



(11) 



such that Op and ep are the initial values at the begin- 
ning of the runaway stage. Therefore, in this case, the 
prospects for planetary ejection are independent of both 
Qfi and the intermediate stages of runaway evolution, as 
long as Qi is high enough to ensure 'I' > 1. Further, as 
in the single stage case, Eq. pip demonstrates that the 
fraction of mass remaining in the star at the moment of 
ejection is equal to (1 -I- ep)/2, independent of the details 
of the intermediate stages. 

If /p = 180°, then the planet's orbit first circular- 
izes before possibly expanding and causing ejection. The 
equations leading to this circularization feature the same 
cancellations, and 



gp (1 + ep) 
- l^-' (1 - ep) 



and 



e. = 1 - (1 - ep) . 



(12) 



(13) 



Hence, prospects for planetary orbit circularization in 
this case are also independent of both ai and the inter- 
mediate stages of runaway evolution, as long as Oi is high 
enough to ensure ^' > 1. The planet becomes circularized 
when the fraction of mass remaining is equal to 1 — ep, 
the same result from Paper I. 

If we consider the mass loss to be instantaneous and 
/p is known then multi-phasic runaway evolution may be 
treated as a consecutive series of impulsive approxima- 
tions (see Section 2.7 of Paper I). We can then create 
relations linking ai, e^ and fi to the same variables from 
all previous runaway stages before ejection. However, be- 
cause the transition between the adiabatic and runaway 
regimes is not sharp (as partially evidenced by the pres- 
ence of k), attempting to link the evolution in both with 
the above methods would likely lead to an unphysical 
result unless the transitional regime is bypassed. 



4.4 Multi-phasic Runaway Regime Evolution 

Investigation of the runaway regime might help assess 
the likelihood of planetary escape when 2> 1. Nonlin- 
ear mass loss in this regime may be treated in a similar 
manner to the adiabatic regime in specific cases. By as- 
suming that mass loss occurs in a series of consecutive 



5 CONCLUSION 

Multiple stars violently interact in ways that a single 
evolving star cannot and the effect on orbiting material is 
hence greater than in the single star case. Conservatively, 
planetary material residing beyond a few hundred AU or- 
biting multiple stars each more massive than the Sun and 
whose minimum pairwise separation is less than 100-Rq 
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is unlikely to remain bound during post-main sequence 
evolution. All Oort cloud analogs in post-main-sequence 
multiple star systems would be disrupted and feature es- 
cape, independent of stellar separations. Planets residing 
at just a few tens of AU from a central concentration of 
stars may be subject to escape in a wide range of multi- 
ple star systems. These systems may provide a significant 
contribution to the free-floating planet population. 

The techniques we utilized in order to obtain these 
results may aid in fu ture studies of individual cases (e.g. 
IVeras fc Wvattll2012l ). We have shown that the prospects 
for planetary escape are determined entirely by the evo- 
lutionary stage with the greatest mass-loss rate relative 
to the remaining system mass. Nonlinear mass-loss rate 
proflles within a phase can be analytically treated by a 
partition into linear segments, where the extent of the 
partition determines the accuracy of the model. For sin- 
gle stars which become white dwarfs, specific constraints 
may be placed on the maximum stellar mass-loss rate 
given assumptions about an orbiting planet. 
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